Spin polarons in triangular antiferromagnets 
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The motion of a single hole in a 2D triangular antiferromagnet is investigated using the t-J 
model. The one-hole states are described by strings of spin deviations around the hole. Using 
projection technique the one-hole spectral function is calculated. For large J/t we find low-lying 
quasiparticle-like bands which are well separated from an incoherent background by a gap of order 
J. However, for small J/t this gap vanishes and the spectrum becomes broad over an energy range 
of several t. The results are compared with SCBA calculations and numerical data. 



Since the discovery of high-temperature superconduc- 
tivity, charge carriers in doped antiferromagnets (AF) 
have been studied intensively. A reliable description of 
the hole motion is important for the understanding of the 
low-energy charge dynamics in the copper-oxide planes of 
the cuprate superconductors. A large number of numeri- 
cal and analytical studies indicate that a single hole in an 
AF spin background has nontrivial properties: The spec- 
tral function consists of a pronounced coherent peak at 
the bottom of the spectrum and a incoherent background 
at larger energies. The coherent peak can be associated 
with the motion of a dressed hole-Le., a hole surrounded 
by spin defects (" spin polaron" ) ETEi. 

Materials with spin arrangements on 2D non-square 
lattices have also been synthetized. For example, ex- 
periments suggest the realization of a triangular .spin-^ 
AF in NaTiO^Q, as well as in surface structuresa such 
as K/Si(lll):B. Delafossite cuprates RCu02+«5, with R a 
rare-earth element, have Cu ions sitting on a triangular 
latticeS. Furthermore, recent results in the context of or- 
ganic superconductors indicate that k-(BEDT-TTF)2X, 
where X is an ion, may be described by a half-filled Hub- 
bard model on an anisotropic triangular latticell3. 

Although most of these materials do not contain a fi- 
nite concentration of holes or electrons away from half- 
filling, it is of fundamental theoretical interest to study 
the hole dynamics in this environment. In this paper 
we adress the motion of a single hole in an otherwise 
half-filled system. The one-hole spectral function for 
this case corresponds directly to the result of an angle- 
resolved photoemission experiment (ARPES) on the un- 
doped compound. 

We assume that a triangular AF doped with holes is 
well described by the t-J model on a triangular lattice, 



H 



(1) 



in standard notation. Note that, in contrast to the square 
lattice, there is no electron-hole symmetry for the trian- 
gular lattice t-J model, i.e., the physics depends on the 
sign of t. For half-filling this model reduces to the tri- 
angular Heisenberg antiferromagnet (THAF) with spins 
ilil. It is by now widely believed thatJts ground state 
possesses magnetic long-range ordert3'Ej which can be 
described as the 120° order being the ground state of the 
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FIG. 1. a) The classical ground state of the triangular AF 
together with the three unit vectors ei of the lattice. The 120° 
spin order defines three sublattices A, B, and C. b) Brillouin 
zone for the triangular system with the path in momentum 
space for the following figures. 

To investigate the hole motion we consider a one- 
particle Green's function describing the creation of a sin- 
gle hole with momentum k at zero temperature: 

G(k,c) = E^ICt^t^IO ( 2 ) 



where z is the complex frequency variable, z — to + irj, 
i] — ► 0. The quantity L denotes the Liouville operator de- 
fined by LA = [H,A]- for arbitrary operators A. \4>o) 
is the ground state of undoped system, i.e., with N elec- 
trons on N lattice sites. 

The one-hole problem in the triangular . latti ce has 
up to now analytically only been studied! 1 ii 15 l using 
self-consistent Born approximation (SCBA). The SCBA 
(without vertex corrections) neglects spiral loops in the 
hole motion (Trugman processes^) since they are formally 
crossing diagrams. However, for the triangular lattice 
these processes are expected to be more important than 
for the square lattice since only three hopping steps are 
necessary for one loop. Here we prefer another analytical 
approach which is based on the picture of the spin-bag 
quasiparticle (QP) or magnetic polaronoQ. The spin de- 
viations surrounding the hole are described by a set of 
path operatorsocl. The correlation-fiinction (g) can be 
evaluated using a cupplant versionEj'EZI of Mori-Zwanzig 
projection techniques. 

In the calculations of this paper the ground state 
l^o^) ofj-thc. undoped THAF is assumed to be long-range 
orderedE3t3. Its description is based on an expansion 
around the classical ground state |(^> c i) of the triangular 
antiferromagnet, see Fig. 1. Translational invariance and 
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rotation invariance in spin space are spontaneously bro- 
ken; we choose the coordinates that the spins in \<p c i) are 
arranged in the x-z plane. 

For the analytical treatment we locally rotate all spins 
to transform the classical ground state into a formally 
ferromagnetic state. The rotation is carried out around 
the y-axis by an angle of Q • with Q = (|tt, 0) at site 
i. The "original" electron operators and the "new" 
operators in the rotated basis are related through: 



Q ■ Ri , . Q Rj 5 

cos — - — d la + sin — - — di-a 
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The resulting Hamiltonian H — Hq + Hi has the follow- 
ing form: 
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are the electron destruction (creation) operators in 
the rotated basis, T are the corresponding spin operators. 
Like in the original t-J model the operators d"' exclude 
double occupancies. The summation symbol (i — > j) 
refers to orientated pairs (ij) where the link i — ► j runs in 
positive direction with respect to one of the unit vectors 
ej of the lattice. The ground state of Hq is now ferromag- 
netic, whereas H± contains the fluctuations. An impor- 
tant difference compared to the square-lattice hole mo- 
tion problem is the existence of a direct hopping term (in 
Ho), i.e., hopping without creation of a background spin- 
defect. This follows from the fact that the spin states on 
adjacent sites in the classical ground state are not exactly 
orthogonal. 

The hole motion process^ will be described in the con- 
cept of path operators auu which create strings of spin 
fluctuations attached to the hole. For the application of 
projection technique we define a set of path operators 
{Ai} which couple to a hole and create local spin defects 
with respect to the classical ground state. The first op- 
erator Aq is the unity operator, the second one Ai moves 
the hole by one lattice spacing creating one spin defect 
and so on. We are interested in calculating dynamical 
correlation functions for the operators {A/Ckc}: 
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G Ia ,j*>(z) = WKA/Ck,,) 



The Green's function (g|) for the physical hole is then 
given by Go^oaf 2 )- Note that it is equivalent to con- 
sider correlation functions for operators {A/dko-}; then 



the physical Green's function (g) has to be calculated 
using relation (^|) . Using cupulants the correlation func- 
tions G can be rewritten astll: 
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( z ) = (<j> cl \n*(A 
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(6) 

The brackets (4> c i\ ■■■ \4>ci) c denote cumulant expectation 
values with \4> c \). The dot • in Eq. indicates that the 
quantity inside (...)' has to be treated as a single entity 
in the cumulant formation. The operator Q transforms 
\4>cl) being the ground state of H into the exact ground 
state \i/jq) of H = Ho + Hi at half-filling. Here is ap- 
proximated with an exponential ansatz which introduces 
spin fluctuations into \4>cl)'- 



= n|</>d) = exp 



(7) 



The operators S v describe the effect of the spin-flip terms 
in Hi. This approach has been shopm to give reason- 
able results for the square lattice AFt3. Here operators 
with up to 4 spin defects with a maximum distance pL 
4 lattice spacings have been employed. Following ref.E£l 
one obtains a non-linear set of equations for the coef- 
ficients a v , = (<f) c \\Sl H Q\<p c \) c , which can be solved 
self-consistently. j_. 

Using Mori-Zwanzig projection tcchniqucEEl one can de- 
rive a set of equations of motion for the dynamical cor- 
relation functions Gi ai j (T /(z). Neglecting the self-energy 
terms it reads: 



/] QKo,Ia{z) Giajgi (z) — XKa,Ja' , 
la 

^Ka,Ja' (z) = z5 KJ 5 a<yl - ^KaXa" Xll" , Ja> • ( 8 ) 

La" 

Xia.ja' and uJia,ja' are the static correlation functions 
and frequency terms, respectively. They are given by the 
following cumulant expressions: 

Xla,Ja> = (0cl|^ t (v4/C k(7 )- t (v4 J c k(7 O'^l0cl) C , 

Ula,Ja> = (0011^(^^)4 (L(Ajc ka ,) )'ft|0 cl ) c . (9) 

These terms describe all dynamic processes within the 
subspace of the Liouville space spanned by the operators 
{Ajdka}. The use of cumulants ensures size-consistency, 
i.e., only spin fluctuations connected with the hole enter 
the final expressions for the one-hole correlation func- 
tion. In the present calculations we have employed up to 
900 projection variables with a maximum paths length 
of 4. The neglection of the self-energy terms leads to 
a discrete set of poles for the Green's functions, in all 
figures we have introduced an artificial linewidth to plot 
the spectra. For details of the calculational procedure 
see e.gJa. 
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FIG. 2. One-hole spectral function for J/\t\ = 0.5 and 
different momenta; left: positive t, right: negative t. The en- 
ergies are measured relative to the energy of a localized hole; 
the path in the Brillouin zone is shown in Fig. 1. 

Now we turn to the discussion of the results. The one- 
hole spectral function Im G(k, to) for both positive and 
negative t and different momenta is shown in Figs. |^ - ^ 
(J/t = 0.5, 2.0, and 10.0). For large J/t we observe pro- 
nounced low-lying QP bands and additional excitations 
with low weight at higher energies which form an inco- 
herent background. With decreasing J/t the gap between 
the QP band and the background excitations decreases. 
For J/t — 0.5 this gap is almost vanished, furthermore a 
sharp QP peak is only present near the bottom of the QP 
band. Values of J/t < 0.2 lead to a incoherent spectrum. 

For negative t the character of the spectral function 
is in principle similar to the t > results, i.e., for large 
J/\t\ (> 1-0) the QP peak is well defined, whereas for 
small J/\t\ the spectra become incoherent. 
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FIG. 3. Same as Fig. 2 but for J/\t\ = 2.0. 



Next we are going to examine the properties of the 
low- lying bands. Fig. || shows the dispersion of the 
lowest pole for positive t and different J/t in comparison 
with data from the literature. For t > the band has 
its energy minimum at momentum (|7r,0) and equiva- 
lent points. For t < and large J/\t\ > 4 the situation 
is similar to the positive t case with the dispersion of the 
QP band reversed, see Fig. ^. However, for intermedi- 
ate J/\t\ and t < the QP dispersion almost vanishes 



(Fig. ||). This can be explained from the composite na- 
ture of the hole motion process in the triangular system. 
The motion consists of direct hopping [amplitude t/2, cf. 
eq. (^[)] and spin-fluctuation-assisted hopping (like in the 
square lattice) with an amplitude being nearly indepen- 
dent of the sign of t (since the main contribution contains 
two hopping steps and one spin- flip process). For nega- 
tive t these two contributions to the dispersion tend to 
cancel each other leading to the very narrow band at 
J/\t\ = 2.0. 
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FIG. 4. Same as Fig. 2 but for J/\t\ = 10.0. 
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FIG. 5. Dispersion of the low-lying QP band for t > 0. 
Left: J/\t\ = 0.5, present calculation (solid) and ED data (12 
sites, crosses) taken from ref.tj. Right: J/\t\ — 10, present 
calculation (solid) together SGBA data (dots) and ED data 
(21 sites, triangles) from ref.cj. The energy zero level has 
been set at the center of mass of the band. 

Another interesting feature is the splitting of the QP 
band into two which is especially visible at large J/\t\. 
It is related to the non-collinear magnetic long-range 
order in the system and can be understood as follows: 
The states |</> c i) and \4>q) have a non-zero helicity, there- 
fore both y directions in spin space (perpendicular to 
the plane defined by the spin directions in |</> c i)) are not 
equivalent. The analysis of the eigenvectors of the dy- 
namic matrix shows that the one-hole eigenmodes 
correspond to states where the missing spin has a definite 
y component. (These are no eigenstates of S** *.) These 
two modes have different energies (for a given k) which 
leads to two distinct bands in the spin-integrated spectral 
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function (g). A spin-resolved photoemission experiment 
(with spins parallel or antiparallel to the helicity direc- 
tion) should observe one or the other of these bands. 

The one-hole dispersion for J/t = 0.5 calculated here 
is in quantitative agreement -with exact diagonalization 
(ED) data on a 12-site clusterEj. The J/t = 10 dispersion 
curve agrees with recent results which where obtained, 
using the SCBA technique and ED on a 21-site clustert£l 
(see Fig. g). To compare results, one has to take into 
account an overall shift of the hole momenta by (|7r,0) 
between clusters with odd and even number of sites. 
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FIG. 6. Comparison of the hole QP bandwidth in triangu- 
lar and square lattice antiferromagnets. Solid/dashed: Tri- 
angular/square lattice, present calculation. Circles: Square 
lattice SCBA results from refB. Inset: J range from to 1. 



Summarizing, we have studied the one-hole motion in 
a triangular AF described by the t-J model. Using the 
spin-polaron concept which describes the hole states in 
terms of local spin defects we find that the one-hole spec- 
tral function shows a QP peak for sufficiently large J/\t\. 
In this regime the picture of a mobile hole dressed by 
spin fluctuations (known from the square lattice) also ap- 
plies to the triangular system. The QP dispersion arises 
both from direct and spin-ffip-assisted hopping processes 
which leads to differences between the t > and t < 
cases. For small J/\t\ the one-hole spectra become inco- 
herent, i.e., the QP weight decreases and the gap between 
the QP band and the background formed by additional 
magnetic excitations vanishes. 

Note that the addition of more than one hole to 
a triangular AF may lead to hole pairing induced by 
spin-wave exchapgeEj which is likely relevant for organic 
superconductorsll3. 

The author thanks K. W. Becker and E. Dagotto 
for useful conversations. Support by the DAAD 
(D/96/34050) and the hospitality of the NHMFL (Talla- 
hassee) are gratefully acknowledged. 



Finally we want to compare the QP bandwidth with 
the one found for the square lattice hole motion. In Fig. 
H we have plotted the bandwidth for both systems de- 
pending _pn the ratio J/t. For the square lattice it is 
knownoB that the bandwidth for small J/t (< 1) is es- 
sentially given by J since hole hopping always creates 
spin defects with respect to the AF background. These 
defects have to be "repaired" by the transverse part of 
the Heisenberg exchange with energy scale J. In contrast, 
in the triangular lattice the bandwidth for small J/t is 
given by at + bJ, whereas for large J/t values it satu- 
rates at 4.5 1. This behavior follows from the existence 
of the direct hopping term with prefactor t/2 in Hq (g) 
which has been discussed above. It leads to a disperison 
proportional to t; the saturation bandwidth is half of the 
bandwidth of the uncorrelated system (9t). Additional 
hole motion processes accompanied by spin fluctuations 
lead to a dispersion proportional to J (at small J/t) as 
in the square lattice. So the physical picture for the hole 
quasiparticle is the following: At small and intermedi- 
ate J/t the hole is surrounded by spin fluctuations, the 
dynamics consists of coherent quasiparticle motion and 
of incoherent processes within the quasiparticle (which 
dominate the spectrum at small J/t). At large J/t spin 
fluctuations are suppressed, but the direct hopping term 
allows the hole to move without creating background spin 
defects. So the pure hole hops with t/2 as a (nearly) free 
fermion. 
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